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Abstract. Context: The shape of dust emission features measured from protoplanetary disks contains information 
about the typical size of the dust particles residing in these disks. A flattened 10 fim silicate feature is often 
interpreted as proof that grain growth has taken place, while a pointy feature is taken as evidence for the pristine 
nature of the dust. 

Aims: We wish to investigate what the effect of dust sedimentation is on the observed 10 /im feature and how 
this may affect the interpretation of the observations. 

Methods: Using a combination of modeling tools, we simulated the sedimentation of a dust grain size distribution 
in an axisymmetric 2-D model of a turbulent protoplanetary disk, and we used a radiative transfer program to 
compute the resulting spectra. 

Results: We find that the sedimentation can turn a flat feature into a pointy one, but only to a limited degree and 
for a very limited set of particle size distributions. If the distribution is too strongly dominated by small grains, 
then the feature is pointy even before sedimentation. If the distribution is too strongly dominated by big grains, 
the sedimentation will not be enough to cause the feature to be pointy. Only if we have a bimodal size distribution, 
i.e. a very small grain population and a bigger grain population, do we find that the transformation from a flat to 
a pointy feature upon dust sedimentation is strong. However, our model shows that, if sedimentation is the sole 
reason for the variety of silicate feature strengths observed in protoplanetary disks, then we would expect to flnd 
a correlation such that disks with weak mid- to far-infrared excess have a stronger 10 ^m silicate feature than 
disks with a strong mid- to far-infrared excess. If this is contrary to what is observed, then this would indicate 
that sedimentation cannot be the main reason for the variety of 10 fim silicate features observed in protoplanetary 
disks. 



Key words, accretion, accretion disks - circumstellar mat- 
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1. Introduction 

The protoplanetary disks surrounding Brown Dwarfs, T 
Tauri stars, and Herbig Ae/Be stars have displayed a 
rich variety of dust emission features in their infrared 
spectra. These features have been studied in great de- 
tail w ith th e Infrared S pace O bservatory (ISO; e.g. Meeus 
at al. I2OOII: Bouwman 120,0 1|) 



with ground-based instru- 
ments (van Boekel et al. i2003i : Honda et al. 2003) and 
with the Spitzer Sp ace Te lescope (e.g. Apai et al. 120051 : 
Kessler-Silacci et al. 120061 : Furlan et al. |2006| ). While the 
wavelength and width of these dust emission features re- 
veals the composition of the dust grains and whether they 
have undergone thermal processing or not, the shape and 
strength of the features betrays the dominant size of the 
emitting grains. Infrared spectroscopy is therefore a pow- 
erful tool for studying the onset of grain growth, which is 



believed to be the very initial step toward planet forma- 
tion. 

A drawback of this kind of analysis is that it only 
probes the very tenuous 'superheated' surface layers of 
the disk. This layer, which is produced by t he irr adiation 
of the disk b y the central star (Calvet et al. Il991 ; Chiang 



tne disK p y tne central star ^tjaivet et ai. iiijijit uniang 
Goldreich|l997|), is very optically thin and has a much 
higher temperature than the disk interior. For these rea- 
sons the dust in this layer typically produce strong emis- 
sion features. Although this surface layer contains only a 
minuscule part of the mass of the disk, it is responsible for 
about hal f of th e infrared luminosity of the disk (Chiang & 
Goldreich ll997r ) and entirely dominates the generation of 
dust emission features. This raises the concern about how 
representative this layer is for the entire disk, and whether 
changes to dust properties observed in this layer are in any 
way related to changes deeper in the disk, where most of 
the mass resides and important processes such as planet 
formation take place. Of particular concern is that dust 
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tends to sediment toward the midplane until an equilib- 
rium between sedimentation a nd ver tical turbulent mixing 
is established (DubruUe et al. Il995l ). Large grains have a 
lower surface-to-mass ratio and are therefore less prone 
to friction with the gas than small grains. They therefore 
settle more deeply into the disk interior before an equilib- 
rium between sedimentation and turbulent mixing is es- 
tablished. Small grains, however, stay afloat much higher 
up in the photosphere of the disk, and since they domi- 
nate the opacity anyway, they maintain the photospheric 
disk height at a higher level. The bigger grains therefore 
settle out of the photosphere, into the deep optically thick 
regions of the disk. They are therefore filtered out of sight. 

The question therefore arises: Could it be that the dust 
in this surface layer is very atypical compared to the bulk 
of the disk deeper in the optically thick parts of the disk? 
If this is the case, what will we learn from studying the in- 
frared spectra from these disks? This paper is not the flrst 
to raise this questi on. It has also been noted by e.g . van 
Boekel et al. |2003|), and Siciha-Aguilar et al. (|2007n have 
invoked the potential effect of big grains sinking out of the 
surface layer in their analysis of Spitzer spectra of a sam - 
ple of pre-main-sequence stars. Kessler-Silacci et al. (l2006l) 
speculate that, by looking at different wavelength ranges, 
one looks deeper into the disk toward regions of larger 
grains. However, no quantitive study of this problem has 
been presented yet. A related problem has, h owever, been 
studied recently by DuUemond et al. (j2007l ). That paper 
showed that dust sedimentation can enhance the infrared 
features from polycyclic aromatic hydrocarbons (PAHs) 
in the disk, because the sedimentation causes the thermal 
dust grains to sink below the photosphere held up by the 
PAHs. 

In the present paper we wish to study the effect of 
sedimentation on the shape of the 10 /im feature of sili- 
cate dust. This solid state infrared emission feature, seen 
in the majority of protoplanetary disks, has often been 
used as a prob e of gr ain size (e.g. van Boekel et al. I2OO5I: 
Schegerer et al. |2006'). A strong and pointy 10 /im feature 
is typical of grains smaller than about 1 /xm in radius, 
while a weak and flattened feature indicates the presence 
of grains between 2 /im and 4/im, clearly larger than the 
averag e grain size in the interstellar medium. Van Boekel 
et al. (|2003h showed that not only are such variations 
in feature shape and strength observed in protoplanetary 
disks, but indeed also the correlation between shape and 
strength expected from the laboratory measurements is 
observed. We are aware that the interpretation of the fea- 
ture shape in terms of grain size is still under debate, as 
porosity and fractal structure of the aggregates may influ- 
ence the ab sorption/emission cross section (Voshchinnikov 
et al. l2006l : Min et al. 2005.). Moreover, Juhasz et al. (A&A 
submitted) has shown that high signal-to-noise is needed 
to interpret these features in the first place. But the fact 
that observations appear to confirm the predicted trend 
gives some hope that the method of interpretation of these 
features is reasonable. 



Many attempts have been made to link the shape and 
strength of the 10 /im feature to the evolution of the disk 
as a whole. But not even the direction of this trend is clear. 
Does the a fiat, broad feature indicate a more developed 
disk by showing that larger grains have been produced in 
the disk? This is the most common interpretation in the 
literature, and it keeps the assumption alive that growth 
from submicron grains to micron-sized grains is a pro- 
cess that happens slowly in the course of a million years. 
However, such a trivial, direct connection seems rather un- 
likely, given the speed at which coagulation can proceed 
and produce large dust grains (DuUemond & Dominik 
l2005h . In view of dust settling, the exact opposite rela- 
tion with the developmental state of the disk is possible 
as well. Assuming a steady-state size distribution and an 
equilibrium state between settling and turbulent mixing, 
the amount of large grains still present in the surface layer 
can be a simple effect of the strength of turbulent vertical 
mixing in the disk. As turbulence and gas densities die 
down toward the end of the disk evolution, large grains 
would tend to disappear from the observer's view and re- 
sult in pointy features for evolved disks. For interpreting 
this type of spectroscopy and linking it to the evolutionary 
state of a disk, it is vital to understand if simple settling 
can indeed filter out large grains efficiently. To establish or 
disprove this possibility is the core purpose of this paper. 

To study the effect of sedimentation on the 10 /im sili- 
cate feature shape, we used the same setup as i n the a bove 
mentioned paper on PAHs (DuUemond et al. l2007l ). We 
set up a disk structure in the same way as described in 
DuUemond & Dominik ( 2004bh . We then assumed some 
grain size distribution throughout the disk and solve the 
time-dependent 1-D vertical dust sedimentation-mixing 
equations at each radial grid point for each grain size. 
Finally we used the Monte Carlo continuum radiative 
transfer code RADMC to produce the infrared spectra, fo- 
cusing on the region between 6 /im and 15 /im. We then 
analyzed these results and figured out if dust sedimenta- 
tion can modify the feature shape. In particular we inves- 
tigated if sedimentation can turn a fiat feature indicative 
of large grains into a pointy feature indicative of small 
grains. 

The structure of this paper is as follows. In Section [5] 
we discuss the opacities we use and how we construct grain 
size distributions. We then discuss the model in Section [3] 
and the results in SectionlH Finally we conclude in Section 

El 



2. Opacities and grain size distributions 

The results of our analysis may depend on the precise 
details of the opacities and grain size distributions used. 
Therefore we discuss them here in some detail. We used 
n and k optical constants from the Jena database for Mg 
Fe Si04. We produced the absorption cross-section using 
a Mie calculation and a specific weight of ^ =3.6 g/cw?. 
We did not include scattering in the radiative transfer 
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Fig. 1. The absorption cross-section in cm^/gram-of-dust 
in the N-band range for the dust we use in our model 
for various spherical grain radii. For grain radii smaller 
than 0.1 iim, the opacity does not change much from the 
opacity of 0.1 fim grains in this wavelength domain. 

calculations. With these Mie calculations we obtained an 
opacity for a grain of a specific radius a. 

To model a grain size distribution we must choose a 
discrete set of grain radii between the minimum radius 
Oniin and the maximum radius flmax of the distribution. 
Throughout this paper we choose Omin — 0.001 /im and 
flmax = 3 /im. In our model we always choose the sam- 
pling points to be logarithmically divided between these 
bounds, such that we have linear spacing in log(a). Figure 
[T] shows the opacities for a limited set of grain sizes, just to 
show the shape change occurring when one goes from small 
(a < 0.3/im) to big a > 2/im grains. It should be noted 
that flat-topped (boxy) features can also be produced by 
certain mixt ures of small crystalline silicate grains (e.g. 



Honda et al. l2003l) . The flatness of the feature is then due 



to the summing of a number of narrow crystalline peaks 
lined up so that, taken together, they look like a broad 
flat-topped feature of an amorphous large grain. High- 
resolution and -sensitivity measurements, however, would 
be able to distinguish these two scenarios by identifying 
the slight bumpy shape of the flat-topped feature in the 
crystalline silicate scenario, where the bumps are at the 
known locations of the peaks of enstatite and forsterite. 

To obtain the total cross section per gram of dust we 
need to speficy the grain size distribution function /(a). 
We define /(a) such that 



Pdust 



m{a) f {a)da 



(1) 



is the total dust density in units of gram/cm^. Here 
m{a) = (47r/3)^a'^ with ^ =3.6 g/cm^. Throughout this 
paper we take a powerlaw distribution: /(a) oc a^. The 
case p = —7/2 repre sents the MRN distribution (Mathis, 
Rumpl & Nordsieck ll977 ). In terms of a size distribution 
in particle mass /(m) oc such that f{m)d'm = f{a)da, 
one gets q = {p ~ 2)/3, so that for the MRN distribution 
we have q = —11/6. However, from now on we shall use 



the /(a) formalism. If we sample this distribution func- 
tion with N discrete a-values, then we can associate a 
dust density to each of these 'bins': 



^dust ~ 



'm{ai)aif{ai)A\og{a)i 



(2) 



where Alog(a)i is the width of each bin in log(a) which 
we assume to be the same for all bins, except for the first 
and the last, which are only half as wide to ensure that 
the distribution goes exactly from log(a,„in) to log(amax)- 
The total radiative cross section ai, in units of 1/cm is 
then 



E 



= 1,N 



The 



total emissivity 



erg s ^ cm ^ Hz ^ ster ^ is then 



units 



(3) 



of 



(4) 



4=1, TV 



where Ti is the temperature of grains with size a^. We 
allow the grains of different sizes to obtain different tem- 
peratures, because we assume that they are not mutually 
thermally coupled, nor do they have strong thermal cou- 
pling with the gas. 

Because each grain size can have a different tempera- 
ture, the a^, is not really a correct representation of the 
expected output spectrum, because the Planck functions 
with which the individual p^ust'^^'*^ multiplied may be 
different for each i. The full radiative transfer calculations 
described later in this paper will take care of this. 



3. Model 

We start our model with making a vertical structure calcu- 
lation of a disk around a Herbig Ae star with T* = lOOOOK, 
i?* = 2.5i?0, M* = 2.5Mq. The disk has an inner ra- 
dius of 0.7 AU, an outer radius of 100 AU, and a sur- 
face density distribution of S(i?) = Eo(^/AU)^^ with 
So — 163g/cm^. We the n use t he procedures described in 
DuUemond & Dominik (l2004al) to create a self-consistent 
disk structure that is in vertical hydrostatic equilibrium. 
This procedure involves solving the radiative transfer 
equation in 2-D axisymmetry using a multi-dimensional 
radiative transfer program, and iterating this with the 
vertical hydrostatic equilibrium equation. This procedure 
then gives pgas{r,z), i.e. the gas density as a function of 
radius and vertical height above the midplane. In this pro- 
cedure we assume that the gas and the dust are well-mixed 
and thermally coupled (i.e. Tgas = Tdust) and we assume 
that the dust consists of small 0.1/xm size silicate grains. 

Once we have this gas density distribution we insert 
the real dust density distribution that we wish to model, 
again with the same dust-to-gas ratio of 1:100. We then 
use the time-dependent dust settling code described in 
detail in DuUemond & Dominik (2004b). This gives the 
density Pi{r, z,t) of each dust component i. For a bimodal 
distribution we have only i = 1,2, while for a continuous 
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distribution we have « = !••• 20, as we model the distribu- 
tion with 20 bins. We take the time t to be 1 Myr, so that 
we are sure that the dust has settled in a settling-mixing 
equilibrium. 

4. Results 

We do the above described experiment for three values 
of the turbulent strength: a = lO^'^, 10^'^, 10"*, and for 
four different continuous size distributions p — —2.5, —2.0, 
—1.0, 0.5. The results are shown in Fig. [2] We find that 
settling can make the feature more pointy, but not dra- 
matically so. This is better seen in Fig. [4]-left, where we 
plot th e well- known strength-shape diagram of van Boekel 
& Min (|2005l ). This diagram shows the strength of the fea- 
ture over the continuum on the x-axis and the normalized 
ratio of the fluxes at 11.3 and 9.8 fim on the y-axis. Small 
grains (pointy features) appear at the lower right of the 
diagram, while big grains (flat features) appear at the top 
right. It should be noted here, however, that crystalliza- 
tion of the grains have an effect similar to grain growth 
by ending up in the top left of the diagram, but we are 
not concerned here about crystallization now. What we do 
expect is that if settling causes a flat feature to be pointy, 
the same model for ever lower values of a would be seen 
as a movement from top left to bottom right. In the fig- 
ure these models are connected by lines. One sees that the 
expected effect is only reasonably strong for p = —1 and 
p = 0.5, but weak or absent for p = —2.5 and p = —2. 

Now we do the same for the bimodal size distribution, 
in which we have only a small grain size (0.001 /im) and a 
big one (3.0/xm). Figure [3] shows that the effect of settling 
is very strong: a flat feature can now really be converted 
into a truly pointy one. For p = —2.8, this happens quite 
suddenly between a = 10~^ and a = lO"'^, while for p — 
—3.25 this happens equally suddenly between a = 10^^ 
and a = 10~^, and for p = —3.5 this happens for any 
value of a. This is also reflected in Fig. |3]-right, where 
one can now clearly see the top-left to bottom-right trend 
predicted when going from flat to pointy. Interestingly, 
such a strong boosting of the small grain population has 
als o bee n seen in the models with PAHs (DuUemond et 
al.|2007|). So it appears that if we have a clearly bimodal 
distribution, then the settling can boost the features of 
the small grains (be they true grains or PAHs), but for a 
continuous distribution this is not as easy. 

In addition to the feature-boosting effect described 
here, dust sedimentation is also known to reduce the 
mid- to far-infrared flux of a disk (Miyake & Na kagawa 
1995t Chiang et al. l200lt DuUemond & Dominik |2004bl : 
To be more precise: the ratio of, 
will decrease for de- 



20061) 



D'Alessio et al 

say, 30 ^m flux over 13 /xm flux 



creasing level of turbulence (and hence increasing level 
of sedimentation). The reason for this is that the disk be- 
comes geometrically flatter as seen in the dust continuum. 
The outer regions of the disk therefore capture less stellar 
radiation and consequently produce weaker mid- to far- 
infrared emission. Depending on the situation, it can even 



happen that the dust in outer disk regions has sunk so 
deeply into the disk that it resides entirely in the shadow 
cast b y the inner disk regions (DuUemond & Dominik 
l2004b[ ). This suppresses the mid- to far-infrared flux even 
more and also makes the disk nearly invisible in scattered 
light. Now, if sedimentation is responsible for the boost- 
ing of small-grain features in some disks, then there should 
be a correlation between the feature strength/shape and 
the mi d- to f ar-infrared flux. In a recent paper, Bouwman 
et al. ( 2008h plot the flux ratio Fao^m/Fia^m versus 10 
fj,m feature-over-continuum for a small sample of T Tauri 
stars. They find a correlation in which the disks with weak 
feature also have a low i^so/jm/^^is/im- They interpret this 
as the effect of grain growth: as the grains grow, the 10 
/im silicate feature becomes weaker and flatter, and be- 
cause the overall opacity of the disk goes down, so does 
the mid- to far-infrared flux compared to the N-band flux. 
However, in the scenario we study in this paper, in which 
different strengths of the 10 /im feature are supposed to be 
due to different strength of turbulence, we would expect 
the opposite correlation: for low turbulence the feature be- 
comes stronger, while the F^o^m/ Fi^^m becomes weaker. 
Indeed this is what comes out of our models, as can be 
seen in Fig. [5] A stronger feature correlates with a smaller 
Fso^m/ Fi3^rn- For the grain size distribution this trend is 
rather weak (but certainly not the opposite trend) , but the 
trend is clear for the bimodal size distribution. If the fea- 
ture strength - SED shape correlation found by Boumwan 
et al. is confirmed on the basis of larger samples, then this 
would serve as clear proof that sedimentation alone can- 
not be held responsible for the strength or weakness of the 
10 /im feature. 



In another recent paper, Watson et al. (|2007r ) found a 
correlation between the crystallinity of the dust as de- 
duced from Spitzer IRS spectra and the SED overall 
shape. Disks that are more crystalline appear to have a 
fiatter geometry. Since there appears to be no physical 
reason why crystalline grains are typically smaller than 
amorphous ones, the sedimentation filtering process we 
discussed in this paper does not predict any correlation 
between these things to happen. The explanation for the 
correlation found by Watson et al. must therefore lie in 
something else. In principle, as mentioned before, a per- 
fect line-up of crystalline peaks in the 9-11 /tm region 
of the spectrum may mimic a flat-topped feature from 
large amorphous sili cates , and may affect the conclusions 
of Bouwman et al. ( 2008h . But it appears that the anal- 
ysis by Bouwman et al. is precise enough to exclude this 
possibility. 

5. Conclusion 

We have modeled the sedimentation of grains of differ- 
ent sizes in a disk, with particular emphasis on the 10 
/im silicate feature shape. We expected that sedimenta- 
tion removes big grains from the surface layers, leaving 
the smaller grains behind, so that the spectrum will then 
be dominated by a pointy 10 /im feature indicative of small 
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Fig. 2. Model results for the continuous size distribution with smallest size 0.001 /im and the biggest size 3 fim. The 
four different panels are for four different values of the powerlaw index p for the grain size. The four different lines in 
each panel are for no settling (top) and for settling-mixing equilibrium with three different levels of turbulence. 
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Fig. 3. As Fig. [21 Model results for the bimodal size distribution with the small grain size 0.001 /im and the big grain 
size 3 ^m. 
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Fig. 4. The feature-strength versus shape diagram of van Boekel & Min (|2005h for both the continuous distribution 
of sizes (left) and the bimodal size distribution (right). Models with the same size distribution are indicated by lines. 
The different symbols are the different levels of turbulence. 
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Fig. 5. The model predictions for the correlation between feature strength (x-axis) and SED shape (y-axis). The SED 
shape is the ratio of the flux (in Jy) at 30 /im and at 13 ^m. Left panel is for the continuous size distribution models. 
Right panel is for the bimodal size distribution models. 



grains, even if the initial size distribution is dominated by 
the flat-topped features characteristic of big (3/im) grains. 
We find that this is indeed the case for bimodal size dis- 
tributions if the smaller grains size is chosen extremely 
small. In our model we have a clearly defined population 
of 3/im size grains and a clearly defined <^ 0.1/im grains. 
The arising effect is reasonably strong because there is a 
large difference in the surface-to-mass ratio between these 
grains, and therefore a strong effect on differential settling. 

On the other hand, our attempt to produce this effect 
with a continuous size distribution have basically failed. 
With a continuous (powerlaw) distribution of grains span- 
ning the range 0.001//m all the way to 3/im, we find 
only a very weak effect, much less than what is observed. 
Therefore, unless the dust size distribution is strongly bi- 



modal distribution with just the right large and small 
grain sizes, this excludes sedimentation alone as the sole 
origin of the different 10 /im silicate feature strengths and 
shapes seen in the spectra of protoplanetary disks. In other 
words, the shape and strength of the silicate feature, which 
originates in the surface layers of the disk, is reasonably 
representative of the dust hidden deep within the disk, 
at least for grains in the size range up to a few micron. 
Changes in the lO/im feature do indicate changes in the 
small grain population in the disk. 

A similar conclusion can be made with respect to the 
influence of the distribution of grains on the overall geom- 
etry of the disk. A pure sedimentation model (i.e. without 
growth) predicts that disks with low 30 /j,m flux over 13 
fim flux ratio should typically have a stronger and pointier 
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10 /im. This is most clearly expected for the bimodal size 
distribution. For the continuous size distribution, this cor- 
relation is expected to be weaker because of the weaker ef- 
fect of sedimentation on the silicate feature shape. But at 
the very least such a model does not predict the opposite 
trend: that weak 30-over-13 flux ratio correlates with a 
weak feature, which is the trend observed by Bouwman et 
al. ()2008). If this trend can be confirmed for a larger sam- 
ple, this would then constitute rather strong proof that 
sedimentation alone cannot be held responsible for the va- 
riety in observed silicate strengths. Dust coagulation, in 
combination of course with sedimentation, must then play 
a vital role in producing the different silicate strengths and 
shapes observed in disks around young stars. 
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